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G\ • Abstract 



The self-gravitating thermal gas (non-relativistic particles of mass m at 
r*** ! temperature T) is exactly equivalent to a field theory with a single scalar 

field </>(x) and exponential self-interaction. We build up perturbation the- 



CN ' ory around a space dependent stationary point 4>o(r) in a finite size domain 



OO 



5 < r < R, (5 << R), which is relevant for astrophysical applications (in- 
. terstellar medium, galaxy distributions). We compute the correlations of 

! the gravitational potential (4>) and of the density and find that they scale; 

the latter scales as r~ 2 . A rich structure emerges in the two-point corre- 
lators from the (j) fluctuations around (j)o(r). The n-point correlators are 
' explicitly computed to the one-loop level. The relevant effective coupling 

turns out to be A = 4ir G m 2 /(T R). The renormalization group equations 
(RGE) for the n-point correlator are derived and the RG flow for the effec- 
tive coupling A(t), r = ln(R/5), explicitly obtained. A novel dependence 
I— i | on r emerges here. A(r) vanishes each time r approaches discrete values 

r = r n = 2-7m/y/7 — 0, n = 0, 1, 2, . . .. Such RG infrared stable behaviour 
[A(t) decreasing with increasing r] is here connected with low density self- 
similar fractal structures fitting one into another. For scales smaller than the 
points T n , ultraviolet unstable behaviour appears which we connect to Jeans' 
unstable behaviour, growing density and fragmentation. Remarkably, we get 
a hierarchy of scales and Jeans lengths following the geometric progression 
R n = Rq e ^n/V7 = Rq [10.749087. . .f. A hierarchy of this type is expected 
for non-spherical geometries, with a rate different from e 2n ^ . 
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I. INTRODUCTION AND RESULTS 



The statistical mechanics of the self-gravitating gas is a relevant physical problem mo- 
tivated by the physics of the cold interstellar medium (ISM) and the large scale galaxy 
distributions [fHH- 

In the grand canonical ensemble, we showed [jXJ] that the self-gravitating gas at temper- 
ature T is exactly equivalent to a field theory of a single scalar field <fi{x) with exponential 
self-interaction through the local euclidean action, 



S[<j>] = ±- fd\ 



~ (V</>(x)) 2 - //V« 



where 



2 /2 „ r /= m Gm 2 
/j, 2 = J- zG VT , T eS = 4tt- 



tt ' T 

Here m stands for the mass of the particles, G for Newton's constant and z is the fugacity 
of the gas. The field equals the gravitational potential up to a multiplicative constant. 
The particle density is given by 

p(r) = --^-V 2 0(r) 

J-eS 

We analyzed this field theory perturbatively and non-perturbatively through the renor- 
malization group approach [0. We showed scaling behaviour (critical) for a continuous 
range of the temperature and of the other physical parameters. We derived in this frame- 
work the scaling relation 

M(R) ~ R dH 

for the mass on a region of size R, and 

Av ~ R q 

for the velocity dispersion where q = |(d# — 1). For the density-density correlations we 
found a power-law behaviour for large distances 

~ |ri-r 2 | M *- 6 . 

The fractal dimension dn is related with the critical exponent v of the correlation length by 

d H = 1/u . 

Mean field theory yields for the scaling exponents v = 1/2, dn — 2 and q — 1/2. Such values 
are compatible with the present ISM observational data: 1.4 < du < 2, 0.3 < q < 0.6 . 

In ref. we developed a field theoretical approach to the galaxy distribution. We 
considered a gas of self-gravitating masses in quasi-thermal equilibrium in the Friedman- 
Robert son- Walker background, We derived the galaxy correlations using renormalization 
group methods. We found that the connected iV-points density correlator C(f\, r 2 , ■ ■ ■ , r^) 
scales as 

' 1 ) 



2 



when t\ » r{, 2 < i < N. The theory is fully predictive and there are no free parameters 
in it. 

Our study of the statistical mechanics of a self-gravitating system indicates that gravity 
provides a dynamical mechanism to produce fractal structure. 

Besides the constant stationary point (f> = — oo studied in ref.( the field theory 
defined by eq . (|1 . 1|) possesses a rotationally invariant and dilatation invariant stationary 
point 

o (r) = ln^l ■ (1-2) 

Non-constant saddle points as 0o( r ) are clearly neccesary in order to describe the physics of 
the self-gravitating gas. 

In the present paper we investigate the perturbation theory around the stationary point 
given by eq.( |1.2j ) in the finite size domain between a very large sphere of radius R and a 
small sphere of radius 6. Such a finite domain 5 < r < R is dictated by the physics of the 
problem. Both in the ISM as well as for the galaxy distributions the relevant fractal domain 
is bounded; beyond R and below S other physics (besides purely gravitational) do intervene. 

The stationary point 0o( r ) may be centered at any arbitrary point x inside the domain 
considered. That is, we must integrate over xo treating it as a collective coordinate 0. 

We explicitly compute the propagator of the field <ft (two points correlator of the gravi- 
tational potential) in the background <po{ r ) to zeroth order [see eq. ( |3.20| )1. It presents a rich 
structure accounting for the fluctuations around </>o(r). 

The n-point correlation functions are explicitly computed at zero momentum to the one- 
loop level. The perturbative series turns out to be an expansion in the effective coupling 

_ g 2 = 4tt G m 2 = T eff 
~ fiR T R R 

Here, 

g 2 = ^T cS ={Mf'^z & ' 2 T ^ (1.3) 

is the dimensionless coupling constant and = dj is the Jeans' length. 

The renormalization group equation (RGE) for n-point correlation functions is derived 
and its coefficients computed to the order A. The RGE gives the variation of the correlation 
functions under a variation of scale. In the present case under variations in R/5. We take 
here into account the explicit dependence on the size of the system through the position of 
the boundaries. 

The new feature in the present model with respect to the standard cases is that the 
propagator has a non-trivial dependence on R/5. As a result, the coefficients in the RGE 
have a novel non-trivial dependence on 

-i R 
r = In — . 

o 

We explicitly obtain the RG flow of the effective coupling A = A(r). Integrating the RGE 
equation from a value r = T{ where A(fj) = Aj is small we find that A(r) decreases until it 
vanishes at the first integer multiple of 2n/\/7. That is at 
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r n = 27rn/v7, n = 0, 1, 2, 



(1.4) 



Just after these points r = r n , the one-loop approximation ceases to be valid since the 
one- loop A(t) becomes negative. 

That is, we can start to run the renormalization group at r = Tj with a small coupling 
Aj = A(Tj) and keep running until r = r n . At this point A(r n ) = 0. In this interval Tj,r„), 
the effective coupling A(r) decreases when the spatial scale r increases as usually happens 
in scalar field theories, which is the case here (infrared stable behaviour) . 

We depict in fig. 1 the running coupling constant A(r) in the intervals (r n _i,r n ) for 
n = 1, 2, 3, 4, 6, 17 as illustrative cases. 

From eq. Ql.4D we see that we have found a hierarchy of scales following the geometric 
progression 

Ro = S, 

R, = R e 2 ^ , 

R n = R e 2nn/V ^ = R [10.749087 . . ] n . (1.5) 

We also expect hierarchies if this type for geometries different from the spherical one. Of 
course, the rate e 2 ^ is expected to change somewhat with the geometry. In addition, 
A(t) grows when r decreases starting from r = n. This suggests that one enters here a 
strong coupling regime. That is, an ultraviolet unstable behaviour that we connect with the 
instability of structures and fragmentation. 

As the effective coupling A(r) grows for decreasing r (ultraviolet unstable behaviour), 
the density fluctuations grow and the Jeans length dj decreases, that is, smaller and smaller 
regions become unstable. We thus get fragmentation of the original mass structure into 
substructures. 

Remarkably, we get a hierarchy of scales and a hierarchy of Jeans lengths following the 
geometric progression 

d Jn = d J0 e 2 ™/ V1 n = 0, 1, 2,... . 

This paper is organized as follows: in section II we summarize the field theory approach 
to the gravitational gas, the properties of the theory under scale transformations and the 
simplest relevant stationary points. In section III we build up the perturbative series around 
the space dependent stationary point 0o( r )- Section IV deals with the renormalization of 
the theory, the derivation of the RGE, the RG flow of the effective coupling A(r), physical 
results and their interpretation. 

II. THE FIELD THEORY APPROACH TO THE GRAVITATIONAL GAS 

As shown in ref. [0,|J, the statistical mechanics of a gravitational gas can be described 
using a continuous field theory description. 
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We consider a gas of nonrelativistic point particles with mass m interacting only through 
their mutual gravitational attraction. We ignore any relativistic or quantum effect. More- 
over, we assume this gas to be in thermal equilibrium at temperature T. We allow the 
number of particles iV to vary through an exchange with the environment. In other words 
we work in the grand canonical ensemble. 

This system can then be described by its grand partition function, namely; 



+00 Z N 



E 

N=0 



N 



j-j- d3 Pl jgU e -(3H N 



1=1 



(2tt)* 



(2.1) 



where (3 = ^ and 



H 



N 2 1 

n = E — - Gm 2 V , 1 , 

^ 2m i<iT1<n hi - ii\ 



(2.2) 



G being Newton's constant and z the fugacity of the gas. 

These N-body expressions can be exactly transformed into a continuum field theoretical 
formula ffl. The number density field is given by 



P(r) 



N 

E^-^O 



and the basic expression for the partition function is : 



Z = Vcj) exp 



T, 



cff 



d 3 * 



;W) 2 -/iV 



where 



(2.3) 



(2.4) 




T, 



cff 



47T- 



Gm 2 



T 



and 



0(x) = 2m^KGf3 ^(x) 



(2.5) 



This expression of the partition function leads to a new insight on the self-gravitating 
gas. Eqs. (|2~l| ) and eq. (|2.4j ) are exactly equivalent. Indeed, the system is exactly described 
by a single scalar field 0(x) with local euclidean action, 



2 ( V< ^( x ))' 



(2.6) 



One characteristic feature of this action is that it is unbounded from below due to the sign 
in front of the exponential self-interaction term. This precisely reflects the presence of an 
attractive force like gravitation. It is then possible that some states evolve into gravitational 
collapse. To avoid the divergencies linked to gravitational collapse we introduce a regulariz- 
ing cutoff at short distances. A simple short distance regularization of the Newtonian force 
for the two-body potential is |j] 



Gm 5 



v„, r 



[l-0( a -r)], 



(2.7) 
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9{x) being the step function. The value of the lower cutoff a depends on the physical problem 
under consideration. 

It is useful to study the mean field equation for the stationary points of the local action. 
It takes the form, 

V 2 o (x) + /iV o(x) = 0. (2.8) 

We give below an equivalence table between the iV-body description and the field theo- 
retical formulation. 



N point particles Continuous scalar field 



U(r) = -Gm £ 7 — > ~m 

i<i<j<N ni-<ij\ m 

N x 



pW^E^r-q,) — -— V 2 0(r) 

Jeflf 



Gm 2 

^— — > T cS 



(2.9) 



A. Behavior under scale transformation 

Let us briefly study the behaviour of the exponential field theory eq.( ^.6[ ) under scale 
transformations. This is a preliminary step to the investigation the renormalization group 
in sec. IV. Scale transformations are defined by 

r^r A = Ar, (2.10) 

where A is a real number. Sometimes we will consider A = 1 + e close to 1. Under T\ the 
field (f) transform as 

0(r)^0 A (r) = 0(Ar) + lnA 2 . (2.11) 
And the classical field equation ( |2.8| ) remains. 



Since (ft is the gravitational potential the addition of the constant In A 2 does not alter the 
physics. 

Our field theory approach can be also generalized to a D dimensional space. The la- 
grangian density £ retain the same structure with D dimensional operators (see M for 
details), 

m = \ (Vd0) 2 - //V, (2.12) 
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where V_d is the D dimensional gradient. 

The scale transformations T\ do not change and the action undergoes the following 
transformation under T\, 

S[cf>] ^ S[<p x ] = X 2 - D S[^]. (2.13) 

The Noether current associated with this transformation is 

J(r) = (r.V0 + 2)V0(r) - rC. (2.14) 

Its divergence in agreement with eq.( |2.13| ), is given by 

VJ = (2 -£>)£. (2.15) 

We did not expect a cancellation of this divergence since the action is not invariant 
under T A except in two space dimensions. One of the important features of this study is 
that, despite the presence of a characteristic length \i in the lagrangian density, the action 
scales under a scale transformation. 



B. Stationary points 

We consider the self-gravitating gas in a large but finite spherical box of radius R. In 
addition, we introduce a short distance cutoff at distance r = 5. This cutoff S ~ a is the 
same as the short- distance cutoff in the gravitational force [eg. ( [2.7|) 1. 

In order to extract information from the functional integral eq.( p.4[ ) we shall search for 
stationary points of the action eq.( [2~6|) . 

Let 0o be a stationary point of the action. This means that 0o is solution of the equation 

V 2 0(r) + //V (r) = (2.16) 

It is interesting to notice that the dilated field (f)\ is also a stationary point. 
The most simple stationary point is |j] 

0o (r) = -oo, (2.17) 

which describes empty space. This is a translationally invariant stationary point. It has 
been studied in ref. |]|. 

Let us now study the non-translationally invariant but rotationally invariant stationary 
point 

o (r) = ln-^ ■ (2-18) 



This solution is also dilatation invariant under the transformation (2.11 



Less symmetrical stationary points certainly exist but their calculation would require 
numerical investigation of the stationary field equation. 
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C. Hydrostatic interpretation of the stationary points 



We recall here that there is a straightforward correspondence between our stationary 
point equation (|2.8|) and the equations for the self-gravitating fluid in hydrostatic equilib- 
rium. In terms of the gravitational potential U(x.) [see eq. (|2.9|) 1, eq. (|2.8| ) takes the form 

V 2 f/(r) = AnGzm i^)^ e~^ (r) . (2.19) 

This corresponds to the Poisson equation for a thermal matter distribution fulfilling the 
hydrostatic equilibrium of an ideal gas. The hydrostatic equilibrium condition is VP(r) = 
— m p(r) VU (r) , where P{r) stands for the pressure, and the equation of state for the ideal 
gas P = Tp yields for the particle density p(r) = p e~^ u ^ , where po is a constant. 
Inserting this relation into the Poisson equation V 2 t/(r) = AttG m p(r) yields eq. (|2.19|) with 



fmT\ 3/2 



The solution 0o given by eq. (|2.18| ) is singular at r = 0. However, such short distance 
singularity is inessential for the long distance properties we are investigating. There exist 
also spherically symmetric regular solutions known as the isothermal sphere H. 



III. PERTURBATION THEORY AROUND A NON CONSTANT STATIONARY 

POINT 

We will build our perturbation theory around the stationary point 0o which is invariant 
under rotations and dilations. That is, 

o (x-x o ) = Inf-aj — fa J • (3.1) 

This stationary point is centered in an arbitrary point xq. 



A. Construction of the perturbation series 

As usual in field theory, we will introduce an external source ^J(x) in order to define 
the correlation functions: 



S[<f>] = ± fee 

J-eS J 



-(V x 0) 2 -/xV 



/ d 3 x ^ J(x) 0(x) . 
J 9 



(3.2) 



J(x) can be interpreted as a test-mass density. Anyway, one sets it equal to zero when 
computing physical quantities. 

Let us now perturb around the stationary point (2.18) 



(x) = o (x - x ) + — x(x - x ) 
v ^ 



(3.3) 
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where 
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cff 



(3.4) 



is the dimensionless coupling constant and x( x ) is? °f course, the fluctuation. The action 
can be then splitted into several terms, 



S[<f>] = S[<f> ] + S 2 [g X ] + S:\gx] - J d 3 x J(x) x (x) , 

where S 2 stands for the quadratic action and Si for higher order couplings, 

1 



d 6 x 



(V X (x)) 2 -/iV«W X 2 (x) 



and 



o g n-2 x „ (x) 



^3 ^~ 2 n\ ■ 



(3.5) 



(3.6) 



(3.7) 



We have a non-polynomial interaction action, hence we shall have vertices with any 
number of legs in the Feynman graphs. 

We compute now the field propagator in the background of the solution 0o( x ) as the 
inverse of the quadratic form in S 2 . 



-V 2 - ;/ 2 >)(|x-x |) 
v x-x () h 1 e 



B. Calculation of the propagator Q in the 4>o background 

We read the inverse propagator Q^ 1 from the quadratic action S 2 , 
£/ _1 ( x - x , y - x ) = £(x - y) 
The propagator Q is defined by the inversion relation: 

J Sz £ _1 (x - x , z - x )£(z - x , y - x ) = 5 3 (x - y). 
We can recast this definition as a partial differential equation, 

-VLx " /.V°(I X - X °D] (?(x - x , y - x ) = 5 3 (x - y). 

Using the explicit form of O and the notation r = x — x ; r' = y — x , we find 

2 



-V z - 



g(r,r')=S(r-r') 



It is convenient to expand the correlation function in partial waves 



m=+l 



1=0 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



where Y^ m {6, <fi) are harmonic spherics and the function gi(r, r') obey the ordinary differential 
equation 



d 2 2d 1(1 + 1) -2 

dr 2 r dr r 2 



gi{r,r') = — S(r - r') . 



(3.13) 



The general solution of the homogeneous differential equation 

d 2 



dr 2 r dr 



2d 1(1 + 1) -2 



yi(r) = 



can be written as 



(3.14) 



A, r ai + B, r 



a, 



where A\ and Bi are arbitrary constants and 



a>i 



'p+i-j 



(3.15) 



As stated in the previous section we consider the gas in a finite volume of radius R with 
a short distance cutoff 5 ~ a. Therefore, we need solutions gi(r,r') in the interval 

5 <r <R, 5 <r' <R. 

The radial propagator gi(r,r') can be then expressed in terms of the solutions of the 
homogeneous equation (|3.14|) as follows: 



gi(r, r') 



1 



r 2 W[yf,yf 



(3.16) 



where 



r< = min(r,r'), r> = max(r,r') 



yf(r) obeys the boundary conditions appropriate for r = 5 and y{"(r) obeys the bound- 
ary conditions appropriate for r = R. W[yf ,yf] stands for the wronskian. Notice that 
r 2 W[yf , yf\ is r-independent. 

The simplest choice of boundary conditions is to set the fluctuations equal to zero both 
at r = 5 and at r = R. That is, 



yf(5) = , y?{R) = 0. 
The radial part of the propagator is now fully determined, it takes the form, 

1 1 



gi(r, r') 



a t sh (a t r) 



sh 











sh 





(3.17) 



(3.18) 



where 



-l R 
t = In — . 

o 



Notice that for I = 0, cto being |eq.( |3.15| )] an imaginary number, the hyperbolic sinus turn 
into trigonometric sines. This is an important point since S-waves produce the leading 
contribution in loop calculations as we shall see in the subsequent sections. 
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9o\r,r 



^sin(^r) Vrr' 



sin 



sin 



2 I P 



(3.19) 



We finally get for the full propagator: 
cj(r.r') 1 ' 2 



47rv / rr 7 [ V7sin(^r 



sm 



sm 



2 I P 



-E 



2/ + 1 



^ a* sh(a ; r) 
where we used the relation 



sh 



ai In 



P 



sh 



En 



oil hi 



2/ + 

47T 



P/(cos7) 



(3.20) 



P z (cos7) . 



and Pz(cos7) stand for Legendre polynomials. 

We can see that the propagator ( |3.25| ) presents a rich structure. It takes into account 
the behaviour of the fluctuations around the background <pQ. Notice that the propagator 
( p.25| ) is symmetric under the exchange of r and r'. 

As one sees from eq.( |3.25| ) for distances such as 5 < r ~ r' < i? the two-point function 
of the gravitational potential scales as ~ 1/ \fr~r' at the tree level. 

Besides if 5 r <C r' <C P the propagator Q(r, r') is dominated by the S-wave. That is, 



nl ,sr,r'»S 1 . (V7 \ . 

e(r ' r) = " sm \~ T sm 



Mi 



sm 



2 I P 



C. Perturbation theory, Generating functionals and Feynman graphs 

We construct now the perturbative series around the stationary point <po- We use stan- 
dard perturbation theory methods ||. We shall give the main steps and the new features 
which appear for our model. 



1. Construction of the general Green functions 

We can write the generating functional of disconnected correlations as follows 

d 3 x 



Z[J] 



AT 



(3.21) 



Z XQ [J] = e J* d3x *>(x-xo).J(x)^ e s 2 [x]-Si[x]+9 f d 3 *. J(x) x(x-xo) (3.22) 

[The name disconnected arises because topologically disconnected Feynman graphs con- 
tibute to the perturbative series of Z[J].] 

Here we used that S[<po] = in three space dimensions. Notice first that the generating 
functional is now a functional of </(.)• F° r J = we recover the previous definition eq. (|2.1| ). 
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As we have already mentioned, 0o( x — x o) is a stationary point for any xo. We have 
therefore to integrate over x . That is, x must be treated as a so called collective coordinate 
0. Such integration over x must be computed exactly since all points x contribute with 
similar weight due to the translational invariance within the large sphere of radius R. In 
order to avoid double counting, the domain of functional integration for the field x( x ) 
must be orthogonal to the linear space of the zero modes. That is, we have to impose the 
constraints 



Jd\ x( x ) d Xi <f) (x - xq) = for i = 1, 2, 3 

M is a normalization factor coming from the Jacobian, 

r 1 " 3 "~ " 3 

7T- jd\ (V0 O : 

_07T J 



(3.23) 



(3.24) 



Since 4>o{x) is rotationally invariant, d Xi 4>o(x), i = 1,2,3 can be expressed as linear 
combinations of the I = 1 modes of the previous section. Namely, yi(r) Yi,±i(9, (ft) and 
yi(r)Yifi(9,(j)). The constraint Q3.23Q implies that we must exclude the I = 1 term in the 
partial wave series. We shall denote from now on, 



Anyr r' 

^ 21 + 1 

f^ 2 o>i sh(a z r) 



>/7sin(^r 



sm 



sh 



sh oil In ( — 



sm 



2 \R . 



Pi(cos'j) 



(3.25) 



It is convenient to split the exponent in eq. (|3.22| ) in two parts, separating the interaction 
part of the action from the gaussian part that can be easily integrated, 

fr)^ e -^fd 3 xxG~ 1 X+Jd 3 xJ(x)x(^.-^o) _ y/fa&g e \j d ' Axdi y J(x)6(x-xo,y-x ) J(y) _ (3.26) 

Using this result and expanding exponentials, one gets after calculation, 

5 



Z X0 [J] = v/det^ef / AJ * £^f-5/ 

n=0 U - \ 



<5J(x) 



~ 1 /l 

X ^Mv2 



fc=0 



d 3 xd\/ J(x) £(x-x ,y-x ) J(y) 



(3.27) 



Z Xo [J] is thus expressed as a power series of J. The coefficients of the ra-th degree term 
is the disconnected n points correlation function of the field (f> (i.e the Green's function). 
That is, 



Q{n) 



1 



Xo 



x l, 



s n z, [J] 



Z xo [0] <fJ(x 1 )...<JJ(x r . 



(3.28) 



j=o 



We will compute and to the first orders in g in terms of Feynman diagrams. 



12 



2. The two points Green function and the Feynman rules 



Let us first compute the sum of "vacuum-vacuum" diagrams (without external sources 

j) z X() [o]. 



2 

Z xo [0] = 1 + f" /^V°^°)[S(x X0 ,x X0 )f 



32 /i 2 



c^ 2 e^ x -M£(xx () ,x X0 )] 3 +0(g t 



(3.29) 



Here we use the notation r xo = r — x . Then, the one and two points correlation functions 
can be written as, 



GS(X)=0O(X XO ) + ^[5] ^(Xxo), 



(3.30) 



Gg(x,y) = 0o(xxo)0o(yxo) + 



^x [0] 



0o(x xo )^ 1) (y X() ) + 0o(yx () )^ 1) (xx o ) + ^ 2) (xx ,yx )] , 

(3.31) 



where, 



V^ 2 



Jd 3 xfi 2 e^o) ^(x xo ,x xo ) S(x xo ,y xo ) 



, /^'7' 2 '' ; - :Xx ' (£(x Xo ,x X() )) 2 £(x Xo ,y Xo ) + ... 

(3.32) 



2 

4 2) (xxo,y Xo ) = £(x xo ,y xo ) + 1-/ rf V 2 e<Mzx(>) ^( x xo> z x ) ^(zxo,Zx ) £(zxo,yx ) + 



(3.33) 



The can be straightforwadly expressed in terms of Feynman diagrams. 



^x [0] 



{1+9+ -I 



9?> = 



{o 



+ 



+...} 
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{ — + +...} 



From these diagrams one can recover the previous integrals using the following Feynman 
rules: 



• For each line: £(x Xo ,y X() ), 

n 2 

The diagram symmetry factor. 



/ \ n— 2 

• For a n-leg vertex: + ( J jd^x/j, 2 e Mxx °\ 



We consider now the connected correlation functions which are simpler than the discon- 
nected ones. 



3. Connected Green functions and vertex functions 
The generating functional of connected correlations is given by 

^ X0 [J]=lnZ X0 [J]. (3.34) 
This is also a series in J whose coefficients are the connected Green functions, 

6 n W X0 [J] 



SJ(x 1 ) . . . 5J(x n ) 



(3.35) 



j=o 



The name connected arises because only connected Feynman diagrams (in the topological 
sense) contribute to the perturbative expansion of the 

The expectation value of the field is given by the one-point function, 

X(x) = -^y. (3.36) 

The generating functional of the vertices (also called 1-particle irreducible) functions is 
defined by the Legendre transformation, 

r X() [x] + ^xo[^]= Jd^xJ (3-37) 
Then, it is easy to derive the inverse relation 

7(x) = (3.38) 



The 1-particle irreducible vertex functions are given as 



rS(xi,...,x n ) 



(3.39) 

J=0 
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Notice that the derivatives are taken at J = 0, that is x — . So, the vertex functions 
are the coefficients of the development of T Xo around a non-zero field. 

The vertex functions are related to the connected Green functions through, 



rg(x) = o, 
r2 ) ) (x,y)=f< ) l~ 1 (x,y), 



(3.40) 



rg(x,y,z) = - Jd% dVi d\ rg>(x, Xl ) rg)( y , yi ) r( 2 o )(z, Zl ) 

Only one-particle irreducible Feynman diagrams contribute to the vertex functions, that 
is, diagrams which cannot be disconnected by a single cut of an internal line. In addition, 
external lines are amputated to these diagrams. 

The knowledge of is crucial since it is the amputating operator. It can be expressed 

as, 



4 2 ) (x,y) = ^ 1 (x,y)-S(x,y), 



(3.41) 



where £(x, y) is the sum of one-particle irreducible diagrams with two amputated external 
lines. In our theory, 

S(x,y) = i/i 2 e^W^(x,x)5(x-y)^ + ^ 2 e ^ x ) (£(x,y)) 2 /i 2 e^ 9 - 

+ 3 2 e0o( x) ( £ (X)X)) 2 J(x _ y) ^ + _ 
Z jji 

or, with diagrams, 



S(x,y) 



+ x + v + 



^ + Oo + •••} 



4- Computation ofT^ andT^ 

We consider in this section I^ 2 ) and T^ 3 ). These are the relevant Green functions for 
renormalization. 



rg = - 2 - <z> - •••} 



r( 3 ) = - 

X() 
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{ ■< + o + <Z> 




where the simple left arrow indicates the inverse propagator. Stumps are meant to remind 
of the amputated external lines. From these diagrams we can write the integrals applying 
our Feynman rules. We have finally to integrate over the collective coordinate x . Moreover, 
choosing to integrate over all but one of the vertex functions variables, yields more compact 
results and keep all the information about the renormalization. We display here the results 
for the first diagrams, where we use the notation 

r = + \n(R/5) , 

As ultraviolet regularization we use a cutoff in the angular momentum sum over I. That 
is, in each propagator we sum up to a maximun L. Since / ~ kr, we use that the momentum 
cutoff is k < 1/8 and the size cutoff is r < R. Therefore, we take 



1<L = R/S 



in each propagator. 



2 

d 9 - J t& |d 3 xo 5(x - y) £(x X() , x X0 ) /i 2 
21 + 1 ch{aiT) 



T V 

l=om 2ai sh(a i r) 



* 21 + 1 

^ a? 

1=0,1^1 "z 



+C 2 t- jdk |rf 3 x £ 2 (x Xo ,y X() ) f e^o) f e *.<y-o) 



9_ 



T 



2 e T 



E 



21 + 1 t {21 + 1) ch(^r) 



+ o E 



2 of sh 2 (^r) 2 , of sh(a z r) 



* 2/ + 1 

=o,/^i "z 



V / = C 5 |rf 3 y |,i 3 x |<i 3 x £(x Xo ,y Xo ) S(y Xo ,z Xo ) S(z Xo ,x Xo ; 



x /i 2 e^^o) /i 2 e 00 ^ ,u 2 e *'( z -o) 



T 



E 



I, (2/ + 1) ch(^r) 3r 



2 e T 



%ti a 3 sh 3 («ir) 2 z= ^ #1 «f sh 2 (a/r) 



+ + £ 



2/ + 1 



.,4„V,2| 



1=0 



3r ^ (2/ + 1) ch(q,r) f * 21 + 1 



2 ,4^! afsh(a,r) 



=o,z^i a i 
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The C{ factors are the symmetry factors of the diagrams. C3 is the factor for the tadpole 
with three amputated external lines instead of two, and C4 is the eye-like diagram with two 
amputated lines on one side and one on the other. Their values are, 

d = C 2 = C 3 = \ , C 4 = \ and C 5 = 1 . 

The only averaged diagrams left to compute, at this order of approximation, are the 
most simple ones, the vertex and the inverse propagator. They yield 



/ d 3 x / d 3 x g 
Jv Jv 



-1, 



-8n(R - 5) , 



(3.42) 



/ rf 3 x / rf 3 x / d 3 y 5(x - yWy - z) -^V ^ = 8n(R - 5) 
Jv Jv Jv ' ^f]x 
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(3.43) 



where V is the volume of the region between the spheres with radii 5 and R. The expansion 
of the zero-momentum vertex functions take the form, 



(r( 2 )> = 



^ J 1 _ -T+lHl 

jv 1 e u^y 8tt 



+- 2 E 



2/ + 1 



+ 1 ch(a ; r) 2 _ 



«=o,/^i 



2af sh(a;r) 



i=fei 2a, 2 sh 2 (a,r) 



1=0,1^=1 a i 



STrfli g I _ T / g 2 \ 1 

J\f y/JIR I 6 V/Ui? J 87T 



E 



I, (2/ + l)ch(a z r) 



+ r 2 v 

3^1 2af sh 2 (^r) 



3(2/ + 1) /t , _. 2 , , _ ^ (2/ + l)ch(a,r) /o , o _ 3 , , 



l + a 2 ) + r £ 



, ,,./, 1 a?sh (a,r) 

^ rti I n/i-r I 

(3 + 3a 2 + af ) 



i=0,< 



E 



2af sh(^r) 



2/ + 1 



+ 6 a l + a l 



(3.44) 



We have now all the quantities to renormalize the coupling constant at one loop order. 
It must be noticed that the perturbative expansion turns out to be in powers of the effective 
coupling 



■ 



(3.45) 



This fact will be very important for the behaviour of the physical quantities. 



5. Zero-momentum value of the T*™) 



Let us investigate the n points vertex functions. We can write at first loop order as 
a sum of rather simple diagrams, 
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r<-> - -{ % + Ci +o< +o +r> + ...i 



n-l vertex n-2 vertex n-2 vertex 



In this summation we should also include all the diagrams deduced from those drawn here 
by topological permutations of the vertex within a given diagram. This would create families 
of diagrams, with each member of a single family containing fcj i-legged vertex. However, 
all members of a given family have the same zero-mementa value (this can be checked on 
the corresponding integrals). 

It is enough for renormalization to compute the diagrams at zero-momenta. Each dia- 
gram drawn above will stand for the whole family and will be assigned the sum Ck s ,...,k n+2 of 
the symmetry coefficients of the members of the family. The averaged zero-momenta value 
of a representative diagram will be noted {{k 3 , . . . , k n+ 2})- We proceed now to compute this 
value. 

Let us call v = Y^ik% the number of vertex in a diagram. Then, 



J^jkj = n + 2v . (3.46) 

3 

We can now recast the diagrams into integrals using Feymann rules and eq. ( |3.46|) , 
({ks,...,k n+2 }) = Cfc 3 ,fc4,...,fc„+2 (~^J fd^od 3 *! ■ • •rf^ti-i^(xi xo ,x 2xo ) • • .£(x VXo ,x lxo ) 



xpVK)...^» 



(3.47) 



The sum Ck 3i k 4 ,...,k n+2 is actually easier to compute than each single symmetry coefficient. It 
can be written as: 



Inserting the explicit propagator in eq.( 3.47 ) and integrating over the angular variables 
yields: 

({k 3 , k 4 , . . . , k n+2 }) = C kzM ^ kn+2 2 v jdn... dr v gi(r u r 2 ) . . . g t (r v , n) 

1=0,1^1 J 

= c k3M _ kn+2 r ±^ a ^ll lT) ( 3 - 49 ) 

where 

v-l 



10 n 



rx rr\ rTv — 1 

(x) = v\ \ dri \ dr 2 . . . / dr v sh 2 (r x + x) TT sh(r fc ) sh(r fc + x) sh 2 (r„) . (3.50) 

Jo Jo Jo j~r 2 
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The coefficients uj v can be directly calculated with formal calculus softwares (we have 
done it for v < 5). In addition, it can be established that they obey the following recurrence 
formula: 



v-l 



Uv[X) = < (-AW) fv[X) + 2^ 

X { i=l L 



(-1) 



i+1 



Ct>„_i UC X" 



where, 



\2^ (ch(x) - 1)' 

{-iTf v+l {x)fr\x)\ 



tanh(f ) ' 



2n 

This recurrence relation allows us to compute the r^ n ^ without performing integrations. 



(3.51) 



However, as far as renormalization is concerned, we need only the leading order of the 
in t. We can first establish that: 



({k 3 , . . . , k n+2 }) =C k3t ... ;kn ^ 



tV * (21 + 1) ch"- 2 (a t r) 
z=t£i «f sh v (a ; r) 



+ r 



v-1 



« £ (21 + 1) ch^Vr) „_ 2 , 



of 1 sh*- i (a,r) 



(3.52) 



Accordingly, only the diagrams with n and n — 1 vertex contribute to the two leading orders 

(ti — 1) ! tl\ 
of T^ n ) in r. Since C n ,o,...,o = ^ an d C n _2,i,o,...,o = ~Ti we find the explicit formula, 



ra-2 



^ )- 



(n-1)! n ^ (2/ + 1) ch"- 2 (a ; r) 
2 T ,4- «r sh B (a,r) 



n! 
+ 4 T 



f (2/ + 1) ch B - 3 (q t r) ,\ n _ 2 l| 



Before we turn to renormalization let us study the density correlation functions. 



D. Density correlations 

The two basic local fields here are the gravitational potential <f)(x) and the particle density 

1 



defined as in eq. (|2.9| ) 



p(x) 



T eff 
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W X) ■ 



The correlation functions of the local density p(x) and their long-range behaviour are par- 
ticularly important. 

Let us introduce a source term in the partition function in order to compute density 
correlations. 

Z p \a\ = JD<P e S[-P]+Jd 3 x P (S)a(x) _ (3 53) 

We then define de density correlation functions as 

'■7 Ul 

(3.54) 



Z p [0] 8a(x 1 ) . . . 5a(x n ) 



(7=0 



The generating functionnal of the connected correlation function of the density field is 
defined as usual as 

W p [a] = In Z p [a]. (3.55) 

We finally define the generating functionnal of the vertex correlation functions of the 
density field: 

T p [p] + W p [a) = J d 3 x p(x) a(x) , (3.56) 

m = ff^ • (3.57) 
oa[x) 

At tree level the two point density correlator takes the form 

< p{x) p(y) > = (^—^j J ^ < V 2 xX (x - So) V 2 x x{y- S ) > 

- 4 />Xq ^( f - f o^- f o) ; ( 3 58 ) 

J\fT eff 2 J \S - x \ 2 \y - x \ 2 

where we used eq. ( |3.11| ) in the course of the calculation. 

Using that (?(x, y) scales as l/J\x\\y\, we see that this correlator scales as l/\x — y\ 2 . 



IV. RENORMALIZATION 

As usual, the renormalized vertex functions are defined as follows: 

r%\ Xl , ...,x n , g R ) = Z n ' 2 rW^, ...,x n ,g). (4.1) 

where Z stands for the wave function renormalization, g is the bare coupling and gn the 
renormalized one. 

The renormalizations are defined such that the two and three points functions take their 
tree level values at zero external momenta. In our case, from eq.( |3.44j ) the renormalization 
prescriptions take the form 
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<rf> 

(Tg>> 



8ttR 

87T.RI 



(1-e-) 



(i-O 



9tl 



(4.2) 



The bare vertex functions at the one- loop level are given by eq. (|3.44j) and can be summarized 
as follows, 



(r< 2 >) 
(r(3)) 



8nR 

aT 

8vri?i 



(i-O 



(1 - O 





C(r)l 




8tt 



1+ ^ 



<? 2 \ D(t) 



fiR J 87T 



(4.3) 



where C(r) and -D(r) can be read off from eqs.( |3.44|) . 

The bare coupling g and the wave function renormalization can be expressed as, 



9r = 9 



1 



1 + 



1 



/ii? y 8-7T 



5 



g 

ui2 



£ 

fiR 

2 \ 2 



(4.4) 



h(r) and -B(r) follow by inserting eqs. (|4.3|) - (|4.4|) into the renormalization prescriptions 
eq.(|4.2|) with the result, 



h{r) = D{t) - -C{t) , B{t) = -C(t) . 

2 

We can thus write Z(r) and up to the first order in 



(4.5) 



Z(t) = 1 + 



fiR) 8vr 



2/ + 1 



— T 



^ (2Z + l)ch(a f r) , 1 2 



1=0 



^ 1 2af sh(a z r) 



+ E 



21 + 1 



' 1 + y) + °(^ 



z=o,z^i 



(4.6) 



0rc = < 1 



S7T 



(21 + 1) ch(a,r) 



r 3 V 

i= w afsh 3 (a z r) 



+ - 2 E 



(2/ + 1) 



b 



0M1 2af sh (ajr 



■(3 



-a, 



+ ^ ^ (2/ + l)ch(^r) 
«=t^i 2af sh(a,T) 



(3 



-a, — 



- E 



=0,1^1 



-^r (4+ 2^-4^ ) + °(^ 



(4.7) 
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A. The renormalization group equations 



1. Derivation of the renormalization group equations 



The renormalization group equations for the vertex functions can be derived as usual, by 
requiring the renormalized vertices to be independent of the scale r except for the explicit 
dependence on the size of the system here arising from the position of the boundary. 

We therefore apply on the bare correlators the differential operator 



dr 2 I dR 85, 



(4.8) 



keeping the bare coupling A defined in eq. fl3.45| ) fixed. 
We have, 

d_ 



rWfx„A. 



T 



Ar (n) . 



where the insertion A arises through the r derivative acting on the bounds of the action: 

1 



A = 



2T, 



eff 



dQ u (R 3 £(Ru) + 5 3 £{5u 



Since we have chosen zero boundary conditions, only the insertion of field derivatives con- 
tribute to the correlators . We can then set 

A = -^r J (RH^xIru) 2 + 5 3 (V X | 5u ) 2 ) dn u . 

Since the boundaries are at finite distances which scale with r, their effect must be taken 
into account through the A insertion. However, on physical grounds the boundary effects 
through the insertion A are expected to be subdominant. In addition, we show in appendix 
B that the insertion has indeed negligeable effects. 

By using the relation between renormalized and bare vertex functions we get, 

3 



9 nf . , d 

a7 + i3(A *W 



n 



- 7 (A 



Hi 



r^(x,,A^,r) = Ar 



(») 

n 



(4.9) 



where we introduced the renormalized effective coupling 

Arc 

The functions (3 and 7 are defined by: 

P(Xn) = 



fiR 



7 (A 



dr 
dlnZ 



■R, 



dr 



(4.10) 
(4.11) 



We find from eq.( |3.53|) that 

7 (A W ) = l + 0(\ n ) ■ 
We now proceed to compute the RG flow of the effective coupling A(r). 
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2. RG flow of the effective coupling A(r) 

The expansion of the effective coupling constant A easily follows from the expansion of 
9, eq.flU, 



A — 

Now, using eq. (|4.10|) , we have: 



1 - Xn^Kr) 



(4.12) 



47r dr 

where h(r) is given by eq. ( |4.4|) and ( |4.7|) 



^ £ («+l)ch(q ; r) +T2 £ (21 + 1) / 3+ 3 a? 
1=0,1^1 a i sh («« r ) z=o,^i 2a i sh («« r ) V 2 

+ 7 £ ( ^.^ ( r ) (« + 3a?-a?)-i £ Hi±i(l6 + 6a?- a f). (4.13) 



4 afsh(a,r) v y 4 /="fei a z 

The renormalization group flow of A is determined by /3(X, r) as usual by 

£=/*(A,r) 
ar 

Notice that contrary to the usual cases, there is a non-trivial and explicit dependence on 
the scale r in h(r), and hence in /3(A, r). We get, 

§ = - J- h'(r) dr. (4.14) 

A z 47T 

Integrating this expression from r = Tj to r, we find for A(r) 

- irx^Ri ■ ^ 

where Aj = A(r^). 

We proceed now to analyse the nontrivial dependence of the effective coupling constant A(r) 
on the scale r. 

Let us define the quantities, 

^21 + 1 



1=2 a i 



We will use the a n for 1 < n < 6: 



/2Z + 1 \ 

o-i = 2(C(0) -!)+£( — 2j = -1.874677 + 0(e 



1=2 



a 2 = 2r - 0.7871527 + 2 e~ T + O (e' 2r ) 
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a 3 = 1.314771 + 2 e~ T + O (e~ 2r ) 

a 4 = 0.4134319 + e" 2r + O (e~ 3r ) 
2 

a 5 = 0.16731628 + - e~ 3r + O (e~ 4r ) 
a 6 = 0.07402368 + — + C (V 5r ) 



Notice that only cr 2 grows logarithmically with the cutoffs. 
This enables us to reduce h(r) to the following form: 

u ^ T " 8 (V7 \ 3 



sin 3 (^r)7V7 V 2 / 49 sin 2 (^r 



37 r /-s/7 \ /3 3 1 

98v/7 sin (^r) V 2 / V2 5 4 3 4 

3 1 10 
-Aa e --a 4 + -(a 2 -2r) - -, (4.16) 



or with the numerical values, 



sin 3 (^r) 7v/7 V 2 / sin 2 (^r) 
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+ tV^cos (—t) - 0.70r-0.55 + C fre- T ) (4.17) 

sin (^r) V 2 / 98^ V ^ 

The effective coupling A(r) exhibits an interesting behaviour as a function of r. 

First of all, the one-loop approximation is reliable only when < A(r) ~ Aj = A(Tj) or 
smaller. Starting at r = , A(r) decreases until it vanishes at the first integer multiple of 

2n/y/7. That is at r n = 2nn/\/7, n = 0, 1, 2, At these points h(r) becomes singular. 

Namely, each time 



sin 



t^Y^j =0 , r = r n = 2nn/V7 



we have A(r n ) = 0. Notice that A(r) becomes negative in the intervals 

T n < t < r n + A„ , 

where A n has a specific value in each interval depending on A« and r« (see fig. 1). [A n 
is defined by the relation: h(r n + A n ) = hfa) — 47r/Aj]. This shows that the one-loop 
approximation ceases to be valid at r = r n + 0. 

The effective coupling constant A(r) decreases with the scale r in an infinite number of 
disjoint intervals 

r„_! + A n _! < r < r n . (4.18) 
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That is, we can start to run the renormalization group at r = Tj, Tj chosen such as 
Tn-i < t,; < r n , with a small coupling Aj = A(rj) and keep running until r = r n . At this 
point A(r n ) = 0. In these intervals Tj,T n ), the effective coupling A(r) decreases when the 
space scale r increases as usually happens in scalar field theories, which is the case here 
(infrared stable behaviour). 

We depict in fig. 1 the running coupling constant A(r) in the intervals given by eq.( |4.18|) 
for n = 1, 2, 3, 4, 6 and 17 as illustrative cases. 

We see that eq.( |4.18| ) provides a hierarchy of intervals where the behaviour of the 
perturbative effective coupling [given by eq.( [4.15| )] is consistent. This hierarchy is formed 
by scales following the geometric progression 

Ro = S, 

R, = R e 2 ^ , 



R n = R e 2nn/V7 = R [10.749087 . . ] n . (4.19) 

Outside the intervals in r defined by eq . ( |4. 18| ) , the effective coupling A(r) becomes neg- 
ative, or very large or both, showing that the one-loop approximation is no more reliable. 



In addition, the growth of A(r) according to eq. ( f4.15|) when r decreases below Tj sug- 
gests that one enters here a strong coupling regime. That is, we find an ultraviolet unstable 
behaviour. We now connect this behaviour with the instability of structures and fragmenta- 
tion. As already established, within a given interval, the coupling A(r) grows with decreasing 
r, that is with a decreasing size of the system (ultraviolet unstable behaviour). In this more 
and more strongly coupled regime, the density fluctuations are ever larger. This means that 
the Jeans length dj decreases, the instability condition (scales larger than dj) becomes easily 
satisfied for smaller and smaller regions, thus leading to the fragmentation of the original 
mass structure into substructures. 

In the opposite regime, A(r) decreases with increasing scale (infrared stable behaviour) 
and the density fluctuations decrease. Jeans' unstabilities thus become less probable as long 
we do not enter the next interval. 

The zeroes of A(r) determine a hierarchy of structures R n = Ro e 2nn ^ described above 
in eq.( |4.19|) . We thus have a self-similar set of structures fitting one into each other. 
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Fig 1: Above is the coupling constant A(r) picture in 6 differents r intervals; [0, ti],[ti, t<i\ 
,[t2, T3],[ts, T4],[t5, tq] and [t\q, Tn]. The initial conditions are set such as Tj = r n — 1 and Aj = 0.01. 
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APPENDIX A: COMPUTATION OF THE MOST SIMPLE FEYNMAN 

DIAGRAMS 



In section 3.3.5 we have given the values of the 2-points and 3-points averaged vertex functions. 
As those values are given for a one loop order computation, only three topologocally different 
amputed diagrams are needed. We will give some details here of the calculation of those diagrams. 



1. Some useful integrals 

The computation of the diagrams, as we will se later, can be brought down to that of three 
integrals. The calculation holds no difficulties, here are the results: 



r-a I 

/ dxsh(x) sh(x + a) = -(ach(a)— sh(a)) 
Jo 2 

r-a i>X I 

J dx J dy sh 2 (x + a) sh 2 (y) = - [a 2 + ach(a) sh(a) — 2sh 2 (a 

a rx-\-a rx 

dx sh(x + a) sh(x) / dysh 2 (y) / dzsh 2 (z) = 
Jo Jo 



Tq ( " Y ch(a) " T sh(a) " \ sh2(a) ch(a) + \ sh3(a) 



Let us see now how these integrals appear in the diagrams. 



2. One loop diagrams calculation 



First is the so-called tadpole diagram, which appear in the 4> 4 theory for example. We compute 
directly its average value. 



2- / = Ci J d\ Jd 3 x <J(x - y) ^(x - x , x - x ) I? e^( x - x °) 

= [ R dr J2(^ + 1 )9l(r,r) 
J 8 i 



E 



21 + 1 f R dr 



ai sh(a;r) J 5 r 



— sh 



a t In ( - 



sh 



r 

ai In ( - 



Let us make the following change in the variables, x = a\ In (|), 



2/ + 1 

, a'f sh(Q;r) 



E.2 



r—aiT 

I dx sh(x + air) sh(a 
Jo 



21 + 1 ch(a ; r) 



E 



21 + 1 



,4^! 2«i sh( a/ r) 



ar 
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The next diagram, that we call eye-like, appear in the c/> 3 theory. 



^/ - -C 2 J d 3 x J d 3 x £ 2 (x - x ,y - x ) /i 2 e ^«( x - x o) f e <A>(y-x„) 
-2^dri ^r 2 ]T( 2/ + l)ftVi,r 2 ) 



We split the domain of integration into two parts where r\ and r 2 are in a given order, either 
r\ < r 2 or n > r 2 . For symmetry reasons the contribution of those parts to the result are the 
same. Beside we implement the same change in variables as in the tadpole diagram, x\ = a; In (^) 
and x 2 = ai In (^), 



= - 4 E 



l= £~jU 1 otj sh (ctyr) n J5 r 2 af sh c^r 
2/ + 1 



r—aiT rxi 
/ dxi / 
JO JO 



,Z^1 a f sh ( a ' r ) 
^ 2/ + 1 



dx 2 sh (xi + ck/t) sh (x 2 ) 



r + Q;r sh(a;r) ch(a/r) — 2 sh (ajr) 



1 2 4-, 2/ + 1 1 * (2/ + l)ch(a;r) 2/ + 1 

" "2 T afsh^r) " 2 T J;^ ofd^r) + 

Last, the triangle diagram, that appears in the 3-points function only, 

= -C 5 Jd 3 y Jd 3 x Jd 3 x g(x - x ,y - x )£(x - y ,z - x ) 

x 0(x - z ,x- x ) /U 2 e^ o(x - xo) m 2 e ^(y" x «) ^ 2 e ^o(z-xo) 

r-R r-R r-R 

= -8^(2Z + l) / dri / dr 2 / dr 3 gi(n,r 2 ) gi(r 2 ,r 3 ) gi(r 3 ,n) 

, J S J S J 8 



R 



Here again we will split the domain into six parts where the three variable are ordered. Contribu- 
tions of theses parts are the same. And using again the same changes in the variables: 
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48 1 -^p-r R ^r-r-^ 2 

i_7TT!ii of sh 6 (aiT) Js r\ Js r 2 Js 



=o,^i a i sh ( ai 



5 r 3 



a x In 



sh 



ct\ In 



sh 



a; In 



r-2 



sh 2 



«z In 



^ (21 + 1) f- a i T f x l f x 2 
~ 48 2] 6 \ / <&1 / dx 2 / ^3 

iJ^i a ° sn \ a i T ) Jo Jo Jo 



sh 2 (xi + ajr) sh (x 2 ) sh (x 2 + ccjt) sh 2 (x 3 ) 
A trick here is to integrate x 2 between and — aor then the integrations on x\ and X3 disentangle, 



= 48 E 



21 + 1 



fei «? sh 3 (a ? r) 7o 



-air 



dx 2 sh(x2) sh(x2+a;T) 



dxi sh 2 (xi) / d,X3 sh 2 (X3) 
Jo 



21 + 1 



3^3 



OCiT 



ah 2 



±— ch(a,T) - - J — sh(a,r) 



sh 2 (air) ch^r) + ^ sh 3 (ccjt) 



* (2/ + l)ch(a ; r) _ 3 2 y. 



2/ + 1 



^ 1 afsh 6 ( ai T) 2 ^^afsh^ajr) 



3 ^ (21 + 1) ch(a ; r) , ^ 2/ + 1 

"9 T Zv ^,5 Q V,^„^\ + Z^ 



- l=Qm ofsh(a,r) 

It is worth to mention that the calculation of any other diagram, averaged or not, is quite 
longer and require the use of formal caculus programs. 



APPENDIX B: CONTRIBUTION OF BORDER TERMS IN THE 
RENORMALIZATION GROUP EQUATION 



To evaluate the r.h.s of the renormalization group equation, we want to compute here the action 
on the correlation functions of the derivative: 



dr 2 V OR 35 



(Bl) 



According to our Feymann rules, a diagram contributing to the n-points vertex function, de- 
noted d XQ (xi, . . . , x n ), with I loops and v vertex writes as follow; 
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dxo {%i > • • • > x 



zi) . . . S(x n - z n )G(zi, ■)■■■ G(-, z v ), (B2) 



where C is the symetry factor of the diagram and the specific form of the propagators product is 
determined by the topological structure of the diagram. We apply now the r derivation and make 
use of our boundary conditions. We find, 



{*i} 



[ — • • • / — 2^( x i ~ z i) ■ • • $( x n — z n )G(zi, .) . . . jr-G(., ■) ■ ■ ■ G{-, z v ), where the sumation is over 

J Z-y J Z v OT 



fd 3 z v 



the derivative positions in the product. We then write: 



d 

■7^G(x,y) = H(x,y) 



(B3) 



where, 



H(r,r') = - V 21 + 1 ■ { di 

A'K\frr' l sh 2 (a;r' 



ol\ In 



6R 



ch (air) — ch 



a\ In • 



|pKcos 7 ) (B4) 



We compare now G and H for r S> 1 and generic values of r and r' when r is not close to its 
critical values; y2 

We can then write | = + ei) , f = \f^0- + e 2)- That is, x and y are not close to the 

borders. Then, 



gi(x,y) = - 



1 1 21 + 1 



1 + (a; - l)ei - (1 + a/)e 2 + 0(ei, e|, eie 2 ) 



(B5) 



hi(x,y) ~ - 



1 2Z + 1 



-raj 



1 + a ; ei + aie 2 + 0(e 1 ,e 2 , eie 2 ) 



In this second equation is defined by analogy with 5;. So, in this region 

gi(x,y) ^ hi(x,y) 1^0 

and 

9o{x,y) ~ h (x,y) 

Then, 

> H(x,y) 



(B6) 



(B7) 



(B8) 



(B9) 



This means that, in this region (i.e not close to the borders), and for non critical values of r, the 
border terms are sub dominant. 



31 



